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Case study — stationary Navier-Stokes solutions 
depending on 11 CAD parameters 


[Audouze, De Vuyst 2008] 
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> Total 11 CAD parameters 
Each Parameter may vary in a small given interval 
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Subdomain of study for POD 


Design of experiment (DoE) — Examples of FE solutions 


Many kinds of solutions 
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| Parameterized discrete FE solutions 


Parameter space : 
u^ = u(x, 0) 


JUN 0 c 10, 1]? 


space vector parameter 


Q : is it possible to get a cheap estimator of 


0 — u^ (.,0) ? 
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Proper Orthogonal Decomposition analysis 
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Comparing a particular solution and its POD 
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Using only K 


(5 POD modes) ! 
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| Plotting all the POD coefficients 
(DoE made of N=150 simulations) - Scatterplot 


Scatterplot coefficients POD 
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| Basics оп POD approximation 


Proposition. Let = > 0 be a small error criterion. If K = K(e) is chosen 


such that m 
> № 
k=K +1 
(Н) DES < €. Proper 
3 M truncation rank 
k=1 
then 
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| 1. POD-Galerkin approach on Laplace problem 


Snapshots: ih _ ,,0°,h 
PN е i uh e yh), 
POD 


W^ = span (L^, we) c y^ 


Variational formulation of Laplace problem: 


(u^. v^) 


> POD-Galerkin 
ГОХ h 0 h h h 
(a^^ yw ln = (fw) va" sew. 
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| Closed form ROM solution and error estimate 


If the Hj-inner product is used for the correlation matrix, then POD modes 
are orthogonal in Hj and 


K 
ROM (х) = S (^, Wh) WH" (x) 
katy 
= ay (0) 
Standard Galerkin error estimate : 
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ui i |е" _ w|] za Yw” = WP. 


Then under (H)-hypothesis, we get 
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| 2. POD-propagator approach 


* Now perform PCA on RHS rather than solutions 
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POD basis : y^ ҮТ q ^.^ c pam 


e [hen build a dual basis such that 


ksh yh kh h h h 
(ers Jay = (Y E Wu P MEUM 
Laplace problem with ХА" 
(m^ as RHS == mode propagator 
= 
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| POD-propagator approach 


Then the solution of the Laplace problem 


(u^ sv") ya = Шш ED vv" c y^. 


is in closed form 
K 


g^ — y (^. ^h oh 


k=1 


Pa Same error estimate as POD-Galerkin 
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| Comparison 
POD-Galerkin POD- 
propagator 
Closed form yes 
ROM 


Intrusive 
approach 


NB of required N K 
Laplace 
solutions oe, 


FOR 


Important aspect in Engineering 
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3. Weighted residual approach 


Remember the Laplace problem : 
(Уи, Vu") — (f*, v^) = 0 Vu" e V^. 
/ i,h . “ : » : 
Let fw |. be the FE “hat function"basis. 


The residual is defined as : 


Ria", f^) = (Уй, Vw") — (f^ wt), 


E Guess 
Ge Residual 
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| Idea 
K 


ü^ (a) = У ap wr? 
k=1 


Find a = a(0) solution of the minimization problem 


"a _ 2 
min 5|[R (u^ (a). f^) [|o 


Numerical complexity far less than a whole 
Laplace problem FE solution 


Non intrusive approach (is the FE code is able to return a residual from a guess) 
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| Non linear problems 


e POD-Galerkin OK 
e POD-Propagator KO (only for linear operators) 
e Weighted residual approach OK 
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4. Data driven approach (non intrusive) 
Parameter space interpolation 


(A (и), v^), = (Р, v") гов) Ур" = ү” 


Looking for а ROM in the form 
K 


a = V а,(0) ve" 


k=l 
Ideally, &ь(Ө) should be 


ак(Ө) = (uh, Wer), 
GC \ Not known 


= 
CENTRALE 20 


P A R I!I & 


| By a design of computer experiment DoCE 


1. Generate a cloud of parameter N points Ө“ In [0, 1]? (Latin Hypercube 
Sampling, Sobol, etc.) ; 


2. Compute N FE solutions и?” = yon 


3. Compute ai, = a;,(6") = (u^^, yeh) 
4. Then interpolate/approximate the a; in the parameter space. 


Э Requires high-dimensional interpolators/approximators. 


Many candidates : 
Low-order polynomials 
Radial basis functions (RBF) [Wendland 2006] 
Kriging approaches 
Moving least squares (diffuse approximation) 


(Fa SPH-like approximation, etc... 
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Case of RBF approximators 


RBF centers 
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| Numerical experiment 


Keller-Segel PDE system 


(modeling skin patterns for living beings) z i Nonlinear terms 


Ou — Au + ау · (uVc) = Bul КИ 


de - вде = (ce 
и 


The system is known to have a rich variety of solutions 
with respect to the different model coefficients. 
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| Exemple of 1D steady-state solutions with periodical ВС 


2 Parameter 
bifurcation 


Here, only y is varying in an interval. 
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Proper orthogonal decomposition 
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POD-RBF ROM solutions compared to FE ones 


(using K=6) 
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| 5. Partition-of-unity local-POD basis ROM 


e POD approximation is sometimes reproached to 
use a basis only dedicated to a particular regime 
(that means for a small parameter-space region) 


Idea > Use different locally-optimal POD basis dedicated to a 
particular parameter region, then reconstruct the global approximation 
in a smooth fashion. 
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| Partition-of-unity local-POD basis ҺОМ 


Parameter domain (0, 1]? 


Practical construction of 
a partition-of-unity: 


Еу 


т, 
O 


D" (0) = Ф | 


Smooth positive compactly 
supported kernel 
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>_ © (0) ш? > 0, РЭ uum = 1. 
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| Algorithm 


1. First, build a local POD basis attached to the 
partition function nb т: 


Km 
min Le, (0°) | yor — Ci T”) yq ^.^ 
Gem — 22€ ON RR o 0T 


(UEP јо, qe Ka 
Bm Weighted least squares 


Km 
Loca ROM: й = У (u^ (Ө), ^) y^^. 


p^ k=l 
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| Global reconstruction 


e Then from the property: 
M 


one can propose the ns reconstruction: 
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-0,h m ) о (Ө Y "(u^ ), v 5^) y^ 
m-—1 


ll 
meim Possible to use 
MoR n again RBF approximators ! 
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| Looking forward : challenges 


e Multilevel ROM modeling for optimization process 


* A posteriori estimators and DoCE enrichment at worst 
error locations for better approximation 


e ROM for nonsmooth solutions (discontinuous solutions 
of hyperbolic equations for example) 


Extension/use in the uncertainty propagation community 
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